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Boost factors of dark matter annihilation into antiprotons and electrons/positrons 
due to the dumpiness of dark matter distribution are studied in detail in this work, 
taking the Sommerfeld effect into account. It has been thought that the Sommerfeld 
effect, if exists, will be more remarkable in substructures because they are colder than 
the host halo, and may result in a larger boost factor. We give a full calculation 
of the boost factors based on the recent N-body simulations. Three typical cases 
of Sommerfeld effects, the non-resonant, moderately resonant and strongly resonant 
cases are considered. We find that for the non-resonant and moderately resonant 
cases the enhancement effects of substructures due to the Sommerfeld effect are very 
small ( ~ 0(1)) because of the saturation behavior of the Sommerfeld effect. For the 
strongly resonant case the boost factor is typically smaller than ~ 0(10). However, 
it is possible in some very extreme cases that DM distribution is adopted to give 
the maximal annihilation the boost factor can reach up to ^ 1000. The variances of 
the boost factors due to different realizations of substructures distribution are also 
discussed in the work. 

PACS numbers: 95.35.+d, 95.85.Ry, 96.50. S- 



I. INTRODUCTION 

The recent observations of the positron excess by PAMELA [1| and the electron excess 

n n 

by ATIC [2[ and PPB-BETS [3| have stimulated extensive discussions. The proposed expla- 
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nations include both the classical astrophysical processes (e.g. 



and possible new physics 



like dark matter (DM) annihilation or decay (e.g. j^). For the annihilating DM scenario a 
"boost factor" of order 1000 of the annihilation rate, compared to that to produce the cor- 
rect relic DM density thermally at the early Universe, is needed to produce enough electrons 
and positrons to fit the data. Although the newly published Fermi result on the electron 
spectrum {gI does not reproduce the ATIC sharp "bump" at 300 — 800 GeV, however, to 
explain the excess at Fermi by DM annihilation requires a similar boost factor 3]- 

There are several mechanisms suggested to generate the boost factor, including the DM 

y, ll^, nonthermal production of DM the Sommer- 
16| ) and the Breit-Wigner resonance enhancement 17|. The 
boost factor due to DM substructures is a natural expectation as N-body simulation shows 
that there are a large amount of substructures exist in the Milky Way (MW) halo 18|. How- 
ever, the detailed calculation based on the N-body simulation shows that the boost factor 
from DM substructures is generally less than 



substructures or DM mini-spi 
feld effect (SE, [l2 
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19|, hereafter we refer it as Paper I). 



Although a single DM clump which is close enough to the Earth may be possible to give large 



boost factor, it is founc 
distribution model 10 



such a case has very small probability to survive in a realistic DM 



20| . Therefore, additional boost effect like the SE or Breit-Wigner 



resonance effect is still necessary in a DM picture to account for the observational data. 

Therefore it is necessary to extend the previous discussions of Paper I on the boost 
factor due to DM clumps to take into account the SE. Furthermore, such a study is also 
very important for the indirect searches of DM. The SE or Breit-Wigner enhancement is 
generally related with the velocity dispersion of DM particles. Comparing with the MW 



halo, the DM in smaller structures is colder and can give larger annihilation signals 2ll. |22|. 



In addition, if the DM substructures can contribute a proper fraction to the locally observed 
electrons/positrons, the boost factor of the Galactic center (GC) can be suppressed due to 
the tidal destroy of substructures in the inner Galaxy, and can avoid the strong constraints 
from 7-ray and radio emission 
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In Refs. 



21 



26|. 



22l | the SE in dwarf galaxies and DM subhaloes has been investigated. How- 



ever, in their works only the effect on the source luminosity is considered. The propagation 
of the charged particles, especially the most relevant electrons/positrons, is not included. 
In this work, we will study the boost factor of DM substructures on the electron/positrons 
and antiprotons after incorporating the SE. 
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To be clear, in the following part of this paper we will refer to the enhancement effect of 
DM substructures with respect to the smooth component as the "boost factor" (see definition 
in Sec. V). For the SE induced enhancement of the smooth component compared with the 
thermal production cross section we will directly call it as "SE enhancement" . 

The paper is organized as follows. We first describe the SE in the next section. Then we 
discuss the configurations of the DM substructures in Sec. III. The propagation models and 
definition of boost factors are simply addressed in Sec. IV. In Sec. V we discuss the results 
of the boost factors from DM substructures including SE. Finally we give the summary and 
conclusion. 



II. THE SOMMERFELD EFFECT 

In the calculation of the DM annihilation cross section at low kinetic energy regime, 
some non-perturbative effects may arise due to new "long-range" interaction. It is called 
"Sommerfeld effect" which requires to sum all ladder diagrams due to the exchange of some 
light scalars or gauge bosons between two incoming DM particles. This effect could lead 
to significant enhancement of the annihilation cross section at small relative velocity. To 
a good approximation, one can simply multiply a factor S to the tree-level annihilation 
cross section. In the calculation of such factor, we use the simplified quantum mechanical 
method in literatures by solving the radial Schrodinger equation with a Yukawa potential 



V{r) = -^e-"'^'' (l3|,ll4| 



— <(r) + -e-™^'-V^(r) = -m^p^tP{r), (1) 



rriy^ r 



where ip{r) is the reduced two-body wave function, and are the masses of DM and 
the new mediating boson respectively, /? is the velocity of DM in the center-of-mass frame. 
In this work, we only consider the simplest situation in which the new mediator is a light 
scalar or Abelian gauge boson. In addition, we only take into account the incoming DM 
^articles with S wave (calculation for arbitrary Ith partial wave could be found in Refs. 



27 



28|). 



As described in Ref. [l^, there are two equivalent methods to achieve S: (1) solving 
the Schrodinger equation with outgoing boundary condition ip' {r) / 'il){r) — ^ im^(3 as r — ^ oo, 
then S is given by S* = |^/^(oo)p/|?/^(0)p; (2) using boundary condition 'ip{r)/r constant 
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as r ^ 0, 5" is given by 5" = |?/''(0)/fcp, where k is the momentum defined as k = m^(3. 
Generally Eq.([T]) needs to be solved numerically. 



In this work, we follow the method of Ref. 



27| to get S numerically. We substitute x = kr 



and (f){x) = Cil){r)/kr into Eq. ([T]) (C is a constant), then solve it with boundary condition 
as 0(0) = 1. We can normalize the solution at infinity as F{x) = x(j){x) C ■sin(x + 5), and 
we have = F'^{x) + F'^{x+7t/2) when x is large enough. Using the formula 5* = (0) / kl"^ , 
the SE enhancement factor is obtained as 



To achieve a more rea. 
of DM particles 



22 



istic result, we also need to take into account the speed distribution 



291 ]. After simply choosing Maxwell-Boltzmann distribution as an 
acceptable approximation, the averaged SE enhancement is given by 

^ = irr-^ r "^p'p" (-^Vi+p^) s{f3), (2) 

J^2[X) Jo ^ 

where x is defined as a; = my./T, K2{x) denotes the second order modified Bessel function of 



the second type, p' = j3 / ^1-/5^ is the normalized DM momentum. When the DM particles 
are non-relativistic in the MW, Eq.(l2]) can be written as 

s = \/f^f"*.^-p(-|^)sM. (3) 

where v and a are the velocity and velocity dispersion of DM respectively, Vesc is the escape 
velocity in the MW. It is obvious that if the behavior of S are S ^ 1/v and S* ~ 1/ v\ the 
average SE enhancement will have similar form S* ~ l/a and S ^ 1/^^ respectively 22 1. 

Some illustrations of the SE enhancements are shown in Fig. [1] In the upper two panels 
the variations of the pre-averaged enhancement factor 5* with model parameters are shown. 
The results are similar with that given in Refs. 2l|, |22| . From the top-right panel we can see 



that there are resonant structures, especially for small velocity (3. If the parameters are tuned 
to approach the resonant region, very large boost factors can be achieved. In the bottom-left 
panel, we show the average SE enhancement factor 5 as a function of the velocity dispersion, 
for several values of m^: close to or away from the resonance peak. In this calculation, we fix 
the mass of DM particle as = 1 TeV, and the coupling constant in the Yukawa potential 
as a = 1/30. Two series of parameter m^, around the two peaks ~ 20 GeV and ~ 1.3 GeV, 
are shown in this panel. We can clearly see the "temperature" dependence of the SE from 
this figure. It is shown that at the early stage of the evolution of the universe, when the 
DM is hot, the SE enhancement is negligible. When DM particles cool down, S become 
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FIG. 1: Top-left: the SE enhancement factor S" as a function of coupling constant a and ra^jm^, 
for a single velocity /3 = 220 km s~^; top-right: S as a function of m^/m-^ for several velocity j3, in 
which a is fixed to be 1/30; bottom-left: the average SE enhancement factor vs. velocity dispersion 
(in unit of light speed) of DM particles, for a = 1/30 and = 1 TeV; bottom-right: the evolution 
of DM abundance with the cosmic time x = m^/T taking into account the SE corresponding to 
the models in the bottom-left panel. 



larger, and finally reach a platform when the DM particles are cold enough. We also label 
the velocity dispersion of today's MW halo, cr ^ 5 x 10~^ (150 km s^^). It can be inferred 
that if the difference of the SE between the smooth MW halo and the saturation value is 
larger, the boost effect from substructures should be more remarkable. We will go details 
on this point in Sec. V. In the bottom-right panel of Fig. [1] we show the evolution of DM 
abundance with the cosmic time x = ttl^/T for the models with = 10, 19 and 19.7 GeV 
shown in the bottom-left panel. We find that there is almost no difference compared with 
the case without SE. The annihilation cross section without SE, i.e. {(Jv)q ~ 1.2 — 1.5 x 10"^^ 
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30| for different models. 



cm^ is found to give the right rehc density of DM 

In the following, we will employ the models shown in the bottom-left panel of Fig. [1] 
to discuss the effects on DM substructures from SE. These models cover the non-resonant 
{m^ = 10 or 1 GeV), moderately resonant [m^ = 19 or 1.24 GeV) and strongly resonant 
{m^ = 19.7 or 1.32 GeV) Sommerfeld enhanced cases without loss of generality. For other 
parameters the conclusions can be easily translated. Note that in this work we will focus on 
the boost factors from DM substructures. We will not go in details of the particle physical 
model of DM or the comparison of the expected fluxes of and p with the data. However, 
some rough implications from the observational data are adopted, like the mass of DM 
^ 1 TeV according to ATIC or Fermi results. Actually as studied in many related 
papers , a DM model with ^ 1 TeV and a total enhancement factor of about 

several hundred can reproduce the observational data. That is to say most of the choosen 
models as shown in the bottom-left panel of Fig. [T] (except for = 10 GeV) are potential 
ones which can explain the data, depending on the dumpiness boost factors that will be 
discussed below. 



III. DM DISTRIBUTION AND SUBSTRUCTURES 



A. Density profile 



The DM density profile based on N-body simulations can be generally parameterized as 
a scale-invariant form 

^ (r/rs)^[l + (r/r3)"](/3-7)/«' ^ 
where ps and are the scale density and radius respectively, and {a, (3, 7) are the shape 
parameters which can be fitted from simulations. The simulations usually favor a central 
cusp of the density profile, like the NFW profile with (a, P, 7) = (1, 3, 1) [3ll] and the Moore 
profile with (a, /5, 7) = (1.5, 3, 1.5) [32], though the exact slope near the center is still under 
debate. In this work we will employ both the NFW and Moore profiles for discussion. Note 
that the central density for these profiles is divergent. To avoid the singularity, we introduce 
a maximum central density Pmax due to the fact that there should be a balance between the 



annihilating rate and the in-fall rate of DM 



For typical parameter settings we have 



10 



18 



10 M0 kpc (Paper I). Throughout this paper we fix pmax to be 10 Mg 
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kpc~^ except special claims. 

The scale parameters ps and are determined using the virial mass M^^^, and the con 
centration parameter Cvk- The virial radius of a DM halo is defined as 

1/3 

,(47r/3)Ap, 

where A ^ 187r2 + 82x-39x2 with x = VlMiz)-! = - 



(5) 



(valid for ACDM universe, 
34| ) is the overdensity, and ^ ISSMq kpc~^ is the critical density of the universe. The 



concentration parameter Cvir is defined as 



r-2 

d(rV) , 

dr \r=r-2 



(6) 



where r_2 refers to the radius at which 
relates rvir and the density profile parameter as (Paper I) 

' vir 

_ Moore ''"vir(^vir 



0. The concentration parameter 



,NFW 



(7) 



Cvir(M,i,)' ' 0.63c,ir(M,ir) 

Therefore if the Cvir — M^i,. relation is specified, rg is determined using Eq.([7]). Finally we 
normalize the total mass J p{r)dV to Mvir to get the scale density ps. 

Generally the Cvir — Afvir relation can be derived through fitting the observational mass 
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However, the observational 



profiles of gravitational systems like galaxy clusters 
sample is limited in a narrow mass range and prevents us from investigating the case down 
to very low mass haloes. Thus similar as in Paper I we will use the toy model predictions 
on the Cv 
ENSOl 



VI ^ Mvir relation based on simulations. Two concentration models BOl 



37| and 



are adopted for discussion, with fitted polynominal form at 2 = as 

4 



ln(cv 



EC. 



X 



i=0 



In 



Me 



with 



and 



a 



BOl 



{4.34, -0.0384, -3.91 x 10"^ -2.2 x 10"^ -5.5 x 10^^}, 



C 



ENSOl 



{3.14, -0.018, -4.06 X 10~^ 0, 0}. 



(8) 



(9) 



(10) 
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B. Substructures 



The mass and spatial distributions of DM substructures can be parameterized as |39l . |4 
diV,ub(r,M,,b 



dPM(M,,b) dPv(r) 
iVsub X — — X 



(11) 



dl^dMsub dMs,b dV ' 

where A'^^ub is the total number of subhaloes, dPn/dMsub and dPv/dy are the normalized 
mass and spatial distribution probabilities. For the mass function, the N-body simulations 
show a power law distribution 



dPM(M 



sub j 



OC 



sub 



(12) 



This relation is assumed to hold in a wide mass range, from the most massive subhalo in 



the MW, M„ 



10^° Mq, down to the scale of Earth mass M^,„ ~ 10"^ M, 







40|. The 



power 
works 



aw index am is about 2, however, with a scattering from ~ 1.7 to ~ 2.1 in various 
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42| . The two recent highest resolution simulations, Aquarius and Via Lactea 



find the values of am to be 1.9 and 2.0 respectively 



of subhaloes will be sensitive^ 



4J]. For am ^ 2.0 the mass fraction 



is still very uncertain 45 



46 



y dependent on the lower cut of the mass of subhalo, which 
4?! . In this work we adopt the fiducial values of am = 1-9 
and the minimum mass Mmin = 10~^ Mq as benchmark model, while the results of other 
parameters are also given for comparison. 

The spatial distribution of subhaloes is usually found to be anti-biased with respect to 
the DM density distribution, and can be fitted with an cored isothermal function (e.g..j39|) 



dPv(r) 
47rr^dr 



oc 



r 



(13) 



where th ~ 0.14r^^^ is the core radius. 

The normalization of the total number of subhaloes is determined by setting the number 
with mass heavier than 10^ M© is 100 (Paper I). For such a normalization and the mass 
function slope am = 1-9, we find the total number of subhaloes with minimum mass M„,\„ = 



10 ^ Mq is about 4 x 10^^, which is consistent with the one obtained in Ref. 



4Q|. 



Finally, we simply give the overall property of the MW halo. We will adopt a NFW 



profile with total mass Mj 



MW 



10^2 Mq 



481 . The virial radius of the MW halo is about 270 



kpc, and the concentration parameter calculated using BOl model is about 13.6. The local 
density is then calculated to be po('"0 = 8.5 kpc) ~ 0.25 GeV cm~^. Note that a fraction of 
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mass / will be in substructures, so the actual density of the so-called "smooth" component 
is (1 — f)p{r). For the above benchmark configuration of subhaloes, / ^ 0.14 is found. Since 
we mainly focus on the boost factors from subhaloes, this smooth halo model is fixed in the 
following discussion. 



IV. PROPAGATION MODEL AND BOOST FACTORS 

In the Galaxy the transport of charged particles is affected by several processes. The 
scattering off random magnetic fields will lead to spatial and energy diffusions. The stellar 
wind may also blow away the cosmic rays (CRs) from the Galactic plane. In addition, 
interactions of CR particles with the interstellar radiation field (ISRF) and/or the interstellar 
medium (ISM) can result in continuous and catastrophic energy losses. Since the detailed 
processes affect the propagation are species-dependent, we will describe the treatments for 
antiprotons and positrons below respectively. The basic common framework is as follows. For 
the transport processes we take a spatial independent diffusion coefficient D{E) = PDqTZ^ 
(where TZ = pc/Ze is the rigidity) and a constant wind Vc directed outwards along z. CRs 
are confined within a cylinder halo L, i.e. the differential density, dN/dE = n, is bound by 
n{z = ±L, -Rmax) = with -Rmax of the scale of the visible Galaxy^. The free parameters of 
the model are the halo size L of the Galaxy, the normalization of the diffusion coefficient 
Kq and its slope S, and the constant galactic wind Vc- 

The propagation equation of CRs can be generally written as 

- D/\N + V,— + 2hT,oAz)N + —i—N] =q{^,E), (14) 

where Ftot = X]i=H He ^« ^ destruction rate of CRs through interaction with ISM in 

the thin gas disk with half height h ^ 0.1 kpc, dE/dt is the energy loss rate, and g(x, E) is 
the source function. 

Given the propagated fiuxes of CRs, we then define the boost factor as the ratio of the 
sum of the smooth and substructure contributions to the smooth one without substructures. 
The detailed formula of the solutions of the propagation equations and the boost factor are 
presented in the Appendix. For more details please refer to Paper I and references therein. 



^ Note that the actual effective volume of CRs depends on the propagation parameters, see e.g., [49 1 
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V. RESULTS 

For the convenience of comparison, we specify the reference model configuration based 
on the descriptions in previous sections: Mmin = 10~^ M0, a^a = 1-9; the inner profile 
of subhalo is NFW, and the concentration model is BOl; the propagation parameters are 
the median settings (see Appendix A). We also include the discussion about an extreme 
configuration with a^a = 2.0, Moore inner profile and BOl concentration model. These two 
parameter settings are compiled in Table [11 

ctm profile Mmin(-^0) concentration propagation 
reference 1.9 NFW 10"^ BOl MED 

extreme 2.0 Moore 10"^ BOl MED 

TABLE I: The reference and extreme model configurations of DM substructures. 

To get a rough idea about how large the enhancement due to SE is necessary to give a non- 
negligible boost factor, we give the relative fluxes of positrons (left) and antiprotons (right) 
in the absense of SE in Fig. [21 It is shown that the contribution to the charged particle 
fluxes of DM subhaloes for the reference configuration is about two orders of magnitude 
lower than the smooth component. If the inner profile of DM subhalo is as cuspy as Moore 
profile, the resulting contribution from DM subhaloes is still about one order of magnitude 
lower than the smooth one. Thus in the absence of SE case, it is very difficult to generate 
large enough boost factor only from DM dumpiness. More details were discussed in Paper 
I. 

In the following we will discuss the cases including the SE. Three kinds of SE enhanced 
cases, i.e., the non-resonant case (m^ = 10 or 1 GeV of Fig. [1]), the moderately resonant case 
{m^ = 19 or 1.24 GeV) and the strongly resonant case (m^ = 19.7 or 1.32 GeV) respectively, 
are discussed one by one. 

A. Non-resonant case 

In Figs. [3land[Hwe show the boost factors for and p respectively, for the non-resonant 
SE model. For the first three panels discussing the effects of DM distribution or propagation 



11 




1 10 100 1000 1 10 100 1000 
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FIG. 2: The (relative) fluxes of positrons (left) and antiprotons (right) from DM annihilation for 
the smooth component (red solid line), the subhalo component with reference configmation (black 
full circle), and the subhalo component with Moore inner profile (blue empty diamond). 



models we set = 10 GeV. Only in the last panel we show the comparison between various 
values of m^. From Fig. [T]we can see that for = 10 GeV, there should be no additional 
enhancement effect on the substructures from SE because the smooth halo has lay in the 
saturation region. In other words the SE enhancement factor will be the same for both 
the smooth component and the substructures. Therefore the results should be the same as 
discussed in Paper I, where we find the boost effects from DM substructures are negligible 
for DM distributions from N-body simulations. 

For the = 1 GeV case, we find a factor about 2 times larger for the saturation value 
than the smooth MW halo. Therefore it will lead to ~ 2 times larger boost factor than that 
for = 10 GeV since most of the subhaloes lie in the satuation region^. This can be seen 
from the extreme case^ in the last panels of Figs. [3] and HI 

From the Figs. [3] and H] we have more details about the boost factor and its uncertainty. 

• Energy dependence of the boost factor B. 



As firstly pointed out in Ref. [20|, the boost factor from DM substructures is energy 
dependent instead of a constant due to the energy-dependent propagation effects of 
charged CRs. This energy- dependent propagation can be further translated into the 



^ For the most massive subhalo A/max ~ 10^° M©, we have cr w 2 x 10 

^ Note, however, since the boost factor is dominated by the smooth component for the reference configura- 
tion, we do not see remarkable differences between = 10 and 1 GeV for the reference model. 
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FIG. 3: The boost factors of positrons for the non-resonant SE case, for different model config- 
urations: varying the distribution of subhalo population (top-left); varying the inner property of 
subhalo (top-right); varying the propagation parameters (bottom-left); and finally the extreme 
cases with Om = 2.0, Moore inner profile with BOX concentration model (bottom-right). The er- 
rorbars show the variances of the boost factors. For clarity of the plot, we slightly shift the energy 
axis among different models in one panel. 



differences of spatial distributions between the smooth component and the substruc- 
tures. This is because the effective propagation lengths of e"*" and p vary with energy 
jiol . I49I . Specificly, for the propagation length decreases with the increase of energy 
due to faster energy loss of high energy positrons; while for p the case is just contrary. 
Thus we may expect a lower boost factor for low energy positrons since it reflects 
the ratio of substructures to the smooth halo in a larger volume, which includes more 
smooth contribution when closing to the GC. Similarly, for antiprotons we will expect 
a higher boost factor for low energy particles. These proterties can be seen clearly in 
the bottom-right panels of Figs. [3] and HI 
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FIG. 4: The same as Fig. [3] but for antiprotons. 



• Energy dependence of the variance ctb. 

The variance of the boost factor ctb is also related to the effective propagation volume, 
and hence the energy, of CRs. The smaller the propagation volume, the larger the 
statistical uncertainty. Therefore the variance is larger for (smaller for p) at high 
energies. These behaviors are also shown in Figs. [3] and HI 

• Dependence on the subhalo mass function. 

Since we use the number of subhalos above 10^ M0 to normalize the number distribu- 
tion function, the lower mass cut Mmin and the mass function slope am will affect the 
fraction of subhalo contribution to the total signals. If Mmin is smaller, or is larger, 
the mass fraction of substructures / will be larger (see Paper I). However, when the 
boost from subhalos is not significant a large fraction of mass in subhalos may lead to 
suppress of DM annihilation, as shown in the figures. This is because the boost factor 
approximates as (1 — /)^ according to Eq. (1B12p . 
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It is not direct to derive the influence on the variance of the boost factor. From 
Eq. flB13|) we know that ctb (Tgub with CTsub given in Eq. flB7l) . The variance of subhalo 
flux will depend on the total number of subhaloes A'sub, the average and variance of 
the annihilation luminosity of a single clump. As shown in Figs. [3] and H] (and also 
the related discussions in Paper I), ctb does not depend sensitively on the parameters 
Mmin and 

• Dependence on the inner property of subhalo. 

In the top-right panels of Figs. [3] and H] we show the boost factors together with the 
variances when changing the subhalo inner profile and the concentration model. It is 
known that BOl concentration model would give larger annihilation luminosity ^ than 
ENS model, and Moore profile would also give larger flux than NEW profile. Therefore 
we get the largest boost factor for Moore + BOl configuration, while the smallest one 
for NEW + ENS configuration, although both are still negligible (~ 1). As for the 
variance, since the annihilation luminosity ^ for ENS concentration model or Moore 
profile differs from the reference model (BOl, NEW) by nearly a constant factor^, we 
expect o-^/{^)m to be almost unchanged. Then from Eq.( lB7l) we have dsub oc ($sub), 
which is the largest for Moore + BOl model and smallest for NEW + ENS model. 

• Dependence on the propagation model. 

The effects of various propagation parameters given in Table [III are shown in the 
bottom-left panels of Figs. [3] and HI Since the propagation parameters affect both 
the smooth component and substructures, the boost factor will not be sensitive to the 
propagation parameters. However, the propagation parameters will strongly affect the 
effective propagation length of particles (A^^ oc \/D^, and oc D/Vc), and then lead 
to different variances. For the minimum parameter setting the effective propagation 
length is the smallest, so the variance is the largest. 

• The maximal case. 

We now come to the extreme model configuration with all the maximum settings of 
parameters discussed above, i.e., BOl concentration model, Moore inner profile and 



It is weakly dependent on clump mass Mgub, see Fig. 3 of Paper I. 
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mass function slope = 2.0^, to show the maximal boost factor. The median 
propagation model is adopted due to the fact that the boost factor does not depend 
on the propagation parameters sensitively. We find that in this extreme case, the 
maximal boost factor will be less than 2 ~ 3. 

B. Moderately resonant case 

If the smooth halo lies out of the saturation region of SE, we may expect stronger boost 
effects from substructures. For the moderately resonant case with m^ = 19 GeV (or 1.24 
GeV), the average SE enhancement factor for the saturation region is ~ 2 (or 4) times larger 
than that of the smooth MW halo. Similar with the previous discussion of = 1 GeV case, 
we will expect ~ 2—4 times larger boost effect than the case without SE. Detailed calculation 
gives exactly the expected results, as shown in Figs. [5] and [61 The detailed properties and 
the parameter-dependences of B and ctb are similar with the previous subsection. We only 
point out that for the extreme case (BOl + Moore + = 2.0) the boost factor is still less 
than, e.g. ~ 5. 

C. Strongly resonant case 

The results for the strongly resonant case with m^ = 19.7 or 1.32 GeV are shown in Figs. 
[7] and [HI According to Fig. [1] we know that the saturation value of the SE enhancement 
factor is about 20 and 10'^ times larger than that of the smooth MW halo for m^ = 19.7 
and 1.32 GeV respectively. For the reference model configuration, we still find almost no 
boost effect for = 19.7 GeV. While for m^ = 1.32 GeV case the final boost factor from 
substructures can be as large as ~ 10 from the bottom-right panel of Figs. [7] and [HI 

For the extreme case of the combinations of all maximal configurations, the maximal 
boost factor is found to be about ~ 20 and ~ 10'^ for m^ = 19.7 and 1.32 GeV respectively, 
see the bottom-right panel of Figs. [7] and [HI We also notice the variance of the extreme case 
if very small. This is because the main contribution comes from the very light microhalos, 
which have large number density. We can conclude that if the SE is finely tuned to be very 



^ When the flux from subhalocs become to dominate the smooth component, the model with am = 2.0 
tends to give the larger boost factor than the one with = 1.8. 
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FIG. 5: The same as Fig. [3] but for the moderately resonant SE enhancement case with = 19 
or 1.24 GeV. 



close to the resonance peak and at the same time taking very cuspy DM density profile and 
subhalo mass function giving huge number of microhalos, the boost effect from substructures 
can be remarkable. However, we think this case is not favorable. On one hand the fine tuning 
of the SE parameter (e.g., in this work) is needed^"^. On the other hand the very cuspy 
subhalo density profile, like the Moore profile, is not favored by the recent high precision 



simulations 



51|. In some of these simulations the density profile is even shallower than 



the NEW profile [511 • We may expect in this case the boost factor due to DM clumps should 
be negligible no matter the SE exists or not. 



^ If we define the degree of fine tuning as 77 = Am^/m^, it is found 77 < 1% is needed to get observable 
boost factors. 

^ Note also that for a very small value of m^, e.g. m^/m^ < 10^^, the saturation velocity can be of 



the order m^/m-^ [50|, which may lead to a non- negligible boost factor even for the non- resonant case. 



However, this can be regarded as another kind of fine tuning. 
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FIG. 6: The same as Fig. [5] but for antiprotons. 



VI. SUMMARY AND DISCUSSION 

It has been shown that the DM dumpiness in the MW halo do not tend to enhance the 
local observed fluxes of charged antiparticles such as antiprotons and positrons, according to 
the N-body simulations of DM structure formation (Paper I). In this work we re-investigate 
this problem taking into account the additional boost effects in DM subhaloes from the SE. 

We flnd that generally the SE, if exists, has the same enhancement effect on the smooth 
component and the substructures due to the saturation behavior, except for flnely tuning 
the SE parameter to some extreme cases when the saturation velocity is much smaller than 
the velocity dispersion of the smooth halo. Therefore for most cases, the conclusions of the 
boost factors from DM dumpiness are similar with Paper I. For the moderately resonant 
case like = 19 or 1.24 GeV in Fig. [H the SE enhancement factor for substructures with 
respect to the smooth halo is also small. The total boost factors from DM subhaloes are also 
negligible for the general cases of the subhalo model (mass function slope < 2.0, BOl 
concentration model, NEW or Moore inner density proflle) given in simulations, as shown 
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FIG. 7: The same as Fig. [3] but for strongly resonant case = 19.7 or 1.32 GeV. 



in the upper panels of Figs. \5\ and O Even for the combinations of several extreme aspects 
of subhalo configurations, e.g., BOl concentration + Moore inner profile + am = 2.0, the 
boost factor is as large as several. 

For the strongly resonant case, e.g. m,^ = 19.7 or 1.32 GeV, the additional enhancement 
factor for substructures from SE is much larger. Even so, for the reference configuration of 
DM distribution which is favored by simulations, we find the total boost factors will be less 
than ~ 10. The largest boost factors for the very extreme case of subhalo configurations 
(BOl concentration + Moore inner profile + ttm = 2.0) can be about ~ 10 to at most ~ 10'^, 
depending on the SE parameter adoption. 

Finally we simply address the implications of DM indirect searches in the GC. As we have 
shown, in general cases the substructures do not play a significant role for the enhancement 
of DM annihilation in the solar neighborhood, as required by the PAMELA and ATIC data. 
Therefore, if the DM interpretation of the observational data is correct, the enhancement to 
the annihilation rate of local DM seems also hold for the GC if there are no other spatial 
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FIG. 8: The same as Fig. [7] but for antiprotons. 



dependent enhancement mechanisms. Thus the photon emission from the GC will be a 
powerful tool to cross check the self-consistency of the theory and constrain the DM density 



in the GC. Alternatively, since the radio and 7-ray emission from t 



le GC in DM annihilation 
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24 
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26|, it may imply an 



scenario is strongly constrained by the observational data 
appeal for a non-negligible contribution to the local DM annihilation from the dumpiness, 
although the DM substructures alone can not fully account for the needed large cross section. 



APPENDIX A: SOLUTIONS OF THE PROPAGATION EQUATIONS OF 
POSITRONS AND ANTIPROTONS 



1. Antiprotons 



It has been shown that for the propagation of antiprotons neglecting the continuous 
energy losses and reacceleration can provide a good enough approach, especially for energies 
higher than several GeV 52|. We will also adopt this approximation here. Therefore the 
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relevant processes include the diffusion, convection and the catastrophic losses — inelastic 
scattering and annihilation in interactions. The propagation equation is 

dN 

- DAN + V,— + 2hV,,,5{z)N = g(x, E), (Al) 

where Ftot = ^i=H He ^« destruction rate of antiprotons in the thin gas disk with 

half height /i ^ 0.1 kpc 521], g(x, E') is the source function. The propagator for a point 



source located at xg, expressed in cylindrical coordinates (r, z) (symmtric in ^) is 52] 



/ 7 \ OO 

gl{r,z,E) = ^^1^^ X Y^c-'K, (r^MTM) sm{KL) sin[fc„(L - z)], (A2) 

n=0 

where r and z are the radial distance and vertical height of the source, Kq{x) is the modified 
Bessel function of the second type, /c„ = Vc/2D, and /c„ is the solution of the equation 
2knCos{knL) = -{2hVtot/D + 2A;„) sin(A;„L), and = 1 - ^^^(kr.L)cos{kr.L) _ somce 
function g(r, z, 9; E), the local observed flux is 

^l{E) = — x2 dz rdrgl{r,z,E) deq{r, z^O; E). (A3) 

4vr Jo Jo Jo 

2. Positrons 

For positrons the case is some different from antiprotons. The dominant process in the 
propagation of positrons is energy loss due to synchrotron and inverse Compton scattering for 
energies higher than ~GeV. In this paper we will neglect the convection and reacceleration 
of positrons, which is shown to be of little effect for E > 10 GeV (also the interested energy 



range here) 53|. Then the propagation equation is 



-DAN + ^(^N] =q{^,E), (A4) 



dE V dt 

in which the second term in the left hand side represents the energy losses. The energy 
loss rate of positrons due to synchrotron and inverse Compton scattering in the MW can 



be adopted as dE/dt = -e^/rE, with e = E/1 GeV and te ^ 10^*^ s [SJ]. We directly 



write down the propagator for a point source located at (r, z) from the solar location with 
monochromatic injection energy E^ |19l . |20| 

g^^ir, E ^ Eg) = ^ X ^o(r, z, f), (A5) 
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in which we define a pseudo time f as 



,5-1 J-l 



f = te ' • (A6) 

^©(r, f) is the Green's function for the re-arranged diffusion equation with respect to the 
pseudo time f 



The effect of boundaries along z = ±L appears in only. Following Ref. 20|] we use two 
distinct regimes to approach Q^^: 

• for ( = L^/ADqt ^ 1 (the extension of electron sphere A = ^^ADqT is small) 

n=— oo V u X / 

where Zn = 2Ln + {—l)"'z; 

• otherwise 

^ oo 

g'^'iz, r) = jJ2 [eM-Doklf)UO)Mz) + expi-Dok'^f)<p'^m'^{z)] , (A9) 
where 



L 

n=l 



Mz) = sm[kn{L - \z\)]; kn = (n - 1/2)tt/L, (AlO) 
c(>'^{z) = sm[k'^{L-z)]- k'^ = nn/L. (All) 

For any source function q{r, z, 6; Es) the local observed flux of positrons can be written as 



V 

4:71 



= ^x2 I dz rdr dEsGI^ {r, z,E^ Es) / d^g(r, z, 6; Es). (A12) 

^ Jo J E Jo 



3. Propagation parameters 

The propagation parameters are determined by fltting the local observed CR spectra 



such as the B/C ratio and unstable-stable isotope ratio. In Ref. 55|] the authors showed 
that there are strong degeneracies between parameters to recover the local measured B/C 
ratio. In this work we adopt three typical settings of parameters which are consistent with 
the B/C data, but give very different primary anti-particle fluxes from DM annihilation 



56(1 . These parameters, labelled as "max", "med" and "min" according to the primary 



antiparticle fluxes, are gathered in Table [TTl 
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S Do 


(kpc^ Myr"^) 


L (kpc) Vc (km s-^) 


max 0.46 


0.0765 


15 


5.0 


med 0.70 


0.0112 


4 


12.0 


mill 0.85 


0.0016 


1 


13.5 



TABLE II: Propagation parameters which are compatible with B/C analysis while giving the 
maximal, median and minimal anti-particle DM fluxes. 



APPENDIX B: BOOST FACTORS 

In this section, we calculate the fluxes of charged anti-particles and the corresponding 
boost factors. The source function of antiprotons and positrons from DM annihilation can 
be written as 



^, (crv) dN 9, , _ ^ , 



where is the mass of DM particle which is set to be ~ 1 TeV in light of the recent 
observations of ATIC and Fermi j^, (av) is the thermally averaged velocity weighted 
cross section, and dN/ dE is the yield spectrum of p or per pair annihilation. As discussed 
in Sec. II, the cross section is velocity-dependent if the SE is taken into account. 

Then we can rewrite the fluxes of p and e~^, Eqs.f lASP and flA12p . for the smooth halo as 



<I>,^{E) = £^ X A X j dWsMOir, 2, E)S{a), (B2) 

where A = (crt')o/2m^, S{a) is the average SE enhancement factor, and Q is the pseudo 
propagator (Paper I) which absorbes the source energy spectrum dN/ dE in the propagator 
of Eq.(jA2| or (jAS]) 

g^{r,z,E) = ^xgUr,z,E), (B3) 



oo 



dN 



g'^{r,z,E)= / dEs-—xgi;{r,z,E^Es). (B4) 
Je (^Es 

Note that in Eq. (lB2p . the SE factor should be spatially dependent due to different velocity 
dispersion of DM particles in the MW. Since the charged particles are thought to come from 
places not very far from us, we will use a© ~ 150 km s^^ to represent the average velocity 
dispersion of the smooth DM halo for simplicity. 
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In our discussion we simplify the source spectrum dN/ dE as in Paper I: for p we adopt 
dN/dE = 1 GeV~^, while for e"*" we set dN/dE = 6{E — m^). This adoption makes the 
following discussion independent of detailed particle physics model of DM, while the major 
property about the boost factor is still kept. It will be easy to convolve any model predicted 
DM source spectrum on these results. 

The flux from DM subhaloes is thought to be the sum of the population of DM point 
sources 

1=1 i=l 

where ^.j = J d^x'p^(x')iS'(crj) is the annihilation luminosity of the zth subhalo. For the 
velocity dispersion of subhaloes, we adopt the relation cr = 2.7 x (Mgub/IO^MQ)^/'^ km s~^ 



22 1 , with the normalization fitted from the observational results of dwarf spheroidals in Ref . 
57|. 

Because we do not know exactly the location and mass of any specified subhalo, there 
should be uncertainty of the prediction of flux from subhaloes. The average and relative 
variance of the flux due to different realization of substructure distribution are (Paper I) 

vA 

($sub) = iVsub X — X (Om X {g)v, (B6) 
47r 



cr. 



sub 



2 




(B7) 



where 



{Om = I dM^^^aMsn^)^^. (B8) 

'AUi„ d^sub 

{g)v = I d^^x^^, (B9) 



V 



max 7D 

4=/, diV«'(M.,b)^-(e)l„ (BIO) 



,2 _ / j3,,/^2d-Pv /^i2 



^1 = j^^'^Q'^-{Qfv. (Bii) 

Finally the boost factor, defined as the ratio of the sum of the smooth and substructure 
contributions to the smooth one without substructures, is* 

B = ^^^^^^^ = (1 - mi + $.ub/f s.), (B12) 

sm 



Note the defination of (f>sm is a bit different from Paper I, where the smooth contribution without sub- 
structures is labelled <f>sm- 
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where $° is the smooth contribution without substructures, and f is the mass fraction of 
subhaloes. The variance of the boost factor is 

ae = = (1 - g^. (B13) 
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